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1. Introduction 

Given the importance of hydrodynamic dynamos in accounting for solar 
and stellar variabilities, a large number of studies have been made of 
their dynamical modes of behaviour using a range of dynamo models. 
However, the complexity and the numerical costs of employing the full 
magneto-hydrodynamical partial differential equations (Oilman 1983) 
have resulted in a great deal of effort going into the study of simpler 
related systems, including mean field and truncated models. 

Both the mean field models (cf. Krause and Radler 1980; Branden- 
burg et al. 1989; Brandenburg, Moss and Tuominen 1989; Tavakol et al. 
1995) and truncated systems (Zeldovich, Ruzmaikin and Sokoloff 1983; 
Weiss, Cattaneo and Jones 1984; Feudel, Jansen and Kurths 1993) have 
been shown to be capable of producing a spectrum of dynamical modes 
of behaviour, including equilibrium, periodic and chaotic states. 

Here we shall concentrate on truncated models and recall that there 
have been two approaches in the study of such models: the quantita- 
tive study of the truncated equations resulting from concrete partial 
differential equations (e.g. Schmalz and Stix 1991) and the qualitative 
approach involving the use of normal form theory to construct robust 
minimal low order systems which capture important aspects of the 
dynamo dynamics (Tobias, Weiss and Kirk 1995). 

Given the approximate nature of the truncated models, an impor- 
tant question arises as to the extent to which the results produced 
by such models are an artefact of their details. This is potentially 
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of importance since, on the basis of results from dynamical systems 
theory, structurally stable systems are not everywhere dense in the 
space of dynamical systems (Smale 1966), and therefore small changes 
in the models can potentially produce qualitatively important changes 
in their dynamics (Tavakol and Ellis 1988; Coley and Tavakol 1992; 
Tavakol et al. 1995). 

To partially answer this question, we take a somewhat complemen- 
tary approach to the previous works, by looking at the robustness of 
the truncated models studied by Schmalz and Stix with respect to a 
number of reasonable changes to their components. This work extends 
and reviews our previous work (Covas et al. 1996). 



2. Mean field dynamo models with a dynamic a 

As an example of how the details of truncated models can change 
their behaviour, we take the truncated dynamic a models considered 
in Schmalz and Stix. The starting point of this work is the mean field 
induction equation 

— = V X (v X B + aB - X B), (1) 

in the usual notation with the turbulent magnetic diffusitivity tjt and 
the coefficient a arising from the correlation of small scale (turbulent) 
velocity and magnetic fields (Krause and Radler 1980). 

For the sake of comparison we take an axisymmetrical configuration 
with one spatial dimension x corresponding to a latitude coordinate 
(and measured in terms of the stellar radius R), together with a longi- 
tudinal velocity with a constant radial gradient (the vertical shear ujq). 
The magnetic field is given by 

B^(0,.4^), ,2) 

where Afp is the (^-component (latitudinal) of the magnetic vector 
potential and Bfp the (/)-component of B. These assumptions allow 
equation @) to be split into 



+ aB^, (3) 



dt Q^2 

dB^ _ 7]t d'^B^ ojQ dA^ 
dt ~ m dx'^ R dx ' ^ ' 

Furthermore a is divided into a static (kinematic) and a dynamic (mag- 
netic) part as follows: a = aocosx — aM{t), with its time-dependent 
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part aM{t) satisfying an evolution equation of the form 

^ = 0{aM)+m), (5) 

where 0{a^j) is the damping term and ^(B) the driving term, which 
is a pseudo-scalar and quadratic in the magnetic field (see Covas et 
al. (1996) for details of this and the truncations of equations (^, 
and (^)). We consider the interval < x < vr (which corresponds to 
the full range of latitudes) and take A = B = C = Oatx = and 
X = TT as boundary conditions. We restrict the allowed modes to the 
antisymmetric subset, that is, those modes that satisfy the condition 
B = at x = tt/2. 

Schmalz and Stix then fix the functional form of J- and study the 
resulting modal truncations. Now given the fact that the exact form of 
J-(B) is not known precisely, we shall examine in the next section how 
robust their results are with respect to changes in J^(B). 



3. Robustness with respect to changes in the driving term 

In their studies, Schmalz and Stix take the following form for the 
feedback term 

.F(B) ~ A^B^, (6) 

and look at the various A^— modal truncations of these equations to find 
out what happens to the dynamical behaviour of the resulting systems 
as A'^ is increased. In all of the following discussion, we assume pure 
antisymmetric modes. 

To determine the nature of the dynamics, we shall employ the spec- 
trum of Lyapunov exponents, and distinguish signatures of the types 
(— , —,—,...), (0, —,—,...), and (-I-, 0, — , . . .) as corresponding to equi- 
librium, periodic, and chaotic regimes respectively. 

For the sake of comparison we summarise, schematically, the results 
of the integration of the systems of Schmalz and Stix (see also Covas et 
al. 1996), in Figure 1. Here the largest Lyapunov exponent is depicted 
by a solid line and its negative, zero and positive values indicate equilib- 
rium, periodic and chaotic regimes respectively (the second Lyapunov 
exponent is plotted as a dashed line). 

Since in many astrophysical settings (including that of the Sun) the 
sign of the dynamo number, D, is not known, we also study the effects 
of changing its sign. The results for negative dynamo numbers are given 
in Figure 2. 

An interesting mode of behaviour occurs in the spiky regions of 
Figure 2 which corresponds to the presence of "multiple attractors" (of 
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D=0 



Dynamo number D 



Figure 1. Schematic graph of the typical asymptotic behaviour of the two largest 
Lyapunov exponents for the D > case with a driving term ^(B) ~ A^B^. The 
route to chaos seems to vary as the truncation order increases. 




Figure 2. Schematic graph of the typical asymptotic behaviour of the two largest 
Lyapunov exponents for the D <Q case with a driving term ^(B) ~ A^B^. 
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Figure 3. Fragility in the dynamics with respect to small changes in the dynamo 
number D and the magnitude of the initial vector {A^, B^, qjv/). Crosses and circles 
represent equilibrium and periodic behaviour respectively. 



more than one attractor consisting of equilibrium and periodic states) 
over substantial intervals of D. 

To clarify the consequences of this behaviour, we have plotted in 
Figure 3 the behaviour of the = 4 truncation as a function of 
small changes in the dynamo number and the initial conditions. It 
clearly demonstrates that in these regimes small changes in either D 
or the initial conditions may produce drastic changes in the dynamical 
behaviour of the system. This form of fragility could be of significance 
in producing seemingly intermittent types of behaviour. 

To study the robustness of this system we considered a modified 
feedback term in the form 

H^) = hA^B^ + ha\-R\\ (7) 

where /i = (/iop) ""^ and /2 = {fiop(3)~^ (/5 is the combined (turbulent 
plus ohmic) diffusion of the field, p the density of the medium and 
the magnetic constant). The additional second term on the right hand 
side of (0) is reminiscent of the type appearing in a more physically 
motivated form of J- given by Kleeorin and Ruzmaikin (1982), Zel- 
dovich, Ruzmaikin and Sokoloff (1983) and Kleeorin, Rogachevskii and 
Ruzmaikin (1995) and discussed in Covas et al. (1996). 
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Figure 4- Schematic graph of the typical asymptotic behaviour of the two largest 
Lyapunov exponents for the case with a driving term ^(B) — fiA^B^ — /2Q|Bp. 
The overall behaviour is independent of the sign of D. 



Figure 4 summarises the results for this modified form of J- and as 
can be seen, the inclusion of a|Bp in equation (|^) has drastic effects 
on the dynamics of the system, for both positive and negative dynamo 
numbers. In particular, it strongly suppresses the chaotic behaviour 
found by Schmalz and Stix. 



4. Comparison with other functional forms 

In Covas et al. (1996), we considered the effects of using more physically 
motivated choices for the driving term given by Kleeorin and Ruz- 
maikin (1982), Zeldovich, Ruzmaikin and Sokoloff (1983) and Kleeorin, 
Rogachevskii and Ruzmaikin (1995) in the form 

T{B)=gi3-B + g2a\B\^, (8) 

where gi = —{fiop)^^ and g2 = {nop(3)~^ are physical constants. 

The results of considering the cases with gi ^ 0, g2 = and gi ^ 0, 
g'2 7^ can be summarised as follows. In the former case, the behaviour 
seems to mirror that seen for the case with ^(B) ~ A(j,B^, in the sense 
that for D > there is an asymptotic dominance of "multiple attrac- 
tor" regimes, while for D < this behaviour is replaced by chaotic 
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behaviour. In the latter case, however, where the term proportional to 
a|Bp is included, there is a suppression of both chaotic and "multiple 
attractor" regimes, and the behaviour looks similar to that depicted 
in Figure 4. We also changed the damping operator so that it was of 
the form derived by Kleeorin, Rogachevskii and Ruzmaikin (1995). It 
was found, however, that this change did not produce any qualitative 
changes. 



5. Conclusions 



We have studied the robustness of truncated a^l dynamos including a 
dynamic a equation, with respect to a change in the driving term of the 
a|Bp type. We find that this changes the results of Schmalz and Stix 
drastically by suppressing the possibility of chaotic behaviour. Our re- 
sults presented here and those in Covas et al. (1996) show that changes 
in the driving term have important effects on the dynamical behaviour 
of the resulting systems. Furthermore, the sign of the dynamo number 
also plays an important role, being capable of radically changing the 
behaviour of the system. Typically we find that the behaviour of the 
system for !F(B) ~ A^B^ and D > is similar to the one with 
~ J-B, but with D < and vice- versa. The change of sign causes 
the behaviour to change from typically chaotic to "multiple attractor" 
solutions. These correspond to regimes where equilibrium and periodic 
regimes are present simultaneously. As a result, small changes in D 
or the initial conditions can substantially change the behaviour of the 
system. This type of behaviour can be of importance in producing 
intermittent type behaviour observed in solar variability. 
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